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We review current cosmological constraints on the axion. We describe the ba¬ 
sic mechanisms by which axions are created in the early universe, focussing on 
the standard thermal scenario where the dominant process is through axion ra¬ 
diation by a string network. A dark matter axion in this case would have a 
mass ma ~ 100 peV, with specified large uncertainties. This cosmological lower 
bound leaves a viable window for the axion below the astrophysical upper limit, 
ma, ^ 1-10 meV. We also discuss alternative axion cosmologies which allow a much 
wider, but indefinite, mass range. 


1 Introduction 

Ever since the early 1980’s, the axion has consistently remained one of the 
most popular dark matter candidates. Unlike many other more exotic parti¬ 
cles, th^axion’s existence depends only on a minimal extension of the standard 
modetl’El which also solves one of its key difficulties—the strong CP problem of 
QCD. In the standard thermal scenario for the early universe, cosmic axions 
are created by a variety of mechanisms ranging from ‘quiescent’ production of 
zero momentum axions when their mass ‘switches on’ at t « 10“® secondsH 
through to ‘topological’ production by the violent ra<hative decay of a network 
of axion strings formed at about t = 10“^® secondafl. Originally, the role of 
axion strings was overlooked and calculations of the present axion density just 
from ‘quiescent’ mechanisms suggested a very light mass near TOq 5 /reV, 
that is, if the axion were to solve the dark matter problerrfl. However, it was 
soon recognized that axion strings provide the dominant cosmological contri- 
butionD’p, and subsequent calculationsQ showed that a heavier mass range was 
appropriate rria ~ 100 /reV, though with additional uncertainties given the dif¬ 
ficulty of these nonlinear calculations. Astrophysical constraints on the axion 
mass suggest that it must be less than about 1-10 meV, so the allowable axion 
parameter window is viable, though relatively narrow. 

The present large-scale axion search experimentLl is looking near the orig¬ 
inal mass range, rria ~ 1-10 peV, for a variety of historical and technological 


1 


reasons. Such an axion is by no means excluded in some alternative scenarios 
which we shall also discuss. However, this paper will focus on the ‘standard 
thermal scenario’ which provides strong motivation for an expanded quest in 
search of a slightly heavier axion—an experiment that may become technically 
feasible in the not-too-distant future. 

1.1 The origin of the axion 

The standard model of particle physics, based around the Weinberg-Salam 
model for electroweak interactions and QCD for strong interactions, has one 
significant flaw—the strong CP problem: non-perturbative effects of instantons 
add an extra term to the perturbative Lagrangian and the coefficient of this 
term, denoted 0, governs the level of CP violation in QCD. The absence of 
any such violation in all observed strong interactions imposes a constraint on 
0, the most strmgent being 9 < 10“^° due to the absence of a neutron electric 
dipole moments. Since the value of 9 is effectively arbitrary, one is left with a 
severe fine tuning problem. 

The elegant solution of Peccei and Quinntl is to allow 6 to become a dy¬ 
namical field which relaxes toward the CP conserving value 0 = 0 H, by the 
spontaneous breaking of a C7(l)-symmetry. The resulting particle, known as 
the axionB, has couplings to ordinary matter which are proportional to f~^ 
and acquires a mass ma also proportional to f~^ at the QCD phase transition. 
Initially, it was supposed that the Peccei-Quinn scale was close to the elec¬ 
troweak phase transition /a ^ Tewi but an exhaustive search of accelerator 
data ruled this possibility out, implying /a^lO’^GeV. However, the ensuing 
disappointment was short-lived because there is no phenomenological reason 
why the Peccei-Quinn scale /a could not be much higher, even up to grand- 
unification scales. Thus, the ‘invisible’ axion was born, an extremely light 
particle with almost undetectably weak couplings of order f~^. 

1.2 Astrophysical eonstraints on the axion 

Accelerator limits on the new ‘invisible’ axion were soon superseded by as¬ 
trophysical calculations of cooling rates of large stars, such as red giants. 
Axions, being weakly coupled, escape from the whole volume of a star and, 
in certain parameter ranges, this can exceed the usual heat loss mechanisms 
through convection and other surface eflfects. The oft-quoted red giant con¬ 
straint is /a^lO®GeV (or ma^6meV)E3. More stringent constraints on /a, 

“In fact, 0^0, since there are some Cl^violating weak interactions, and it has been 
shown that 9 > in the standard modejj, although this is highly model dependent. 
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however, were thought to be imposed by neutrino emissions from supernova 
1987a, provic^g the limit, /a^lO^°GeV, which corresponds to a mass bound 
0.6 meVllil. This ‘conventional’ wisdom has been questioned by Raffelt et 
al. who performed a detailed re-examination of these astrophysical calculations. 
They suggest that one cannot impose a meaningful constraint from SN1987a, 
since there is insufficient data, and, moreover, that the constraint from md 
giants is weaker than previously thought, /a ^ 5 x 10®GeV (or rua ^ 10 meV)llj. 
More recent work endeavours to answer some of these criticisms and reasserts a 
somewhat weaker bound from SN1987a, /a^lO^GeV (or ma^6meV), though 
this can Jae an order of magnitude tighter depending on the axion model pa- 
rametersEj. This topic lies outside t^ scope of the present discussion, so more 
general reviews should be consultedEJ. 


1.3 Cosmological axions: ‘quiescent’vs. ‘topological’production 

The weak coupling to ordinary matter and the substantial redshifting between 
the Peccei-Quinn and QCD phase transitions make the axion an ideal cold d^'k 
matter candidate. The earliest estimates of the cosmological axion densit3S0 
assumed ‘quiescent’ production, that is, axions were created through coherent 
oscillations about the minimum of the Peccei-Quinn potential when the axion 
mass ‘switched on’ at T ~ Aqcd- A recent estimate of the relative contribution 
of these zero momentum axions is 




a,h 


0.9/i"^A 


/a \ 

IQi^GeV ) 
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where A ~ 1 ^counts for the uncertainties in the QCD phase transition (dis¬ 
cussed later)EA. Assuming an ill-defined ‘average’ value of the axion field 
m) = 7r^/3, implies a constraint /a^ lO^^GeV, Wa^S/reV when compared 
to the closure density of a flat universe fl = 1. As we have emphasised already, 
this estimate ignores much stronger topological effects arising at the Peccei- 
Quinn phase transition, that is, the inevitable creation of an axion string net¬ 
work. 

If we suppose that we live in a universe that underwent an inflationary 
phase in the early universe, then there are two basic scenarios for cosmological 
axion production. These depend on the relationship between the inflationary 
reheat temperature Tj-eh and the string-forming phase transition temperature 
Tpq - U 

1. If the reheat temperature T^eh is sufficient to restore the Peccei-Quinn 
symmetry Treh > /a, then a network of global or axion strings will form 
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by the Kibble mechanism. This is the usual thermal scenario commonly 
assumed in most of the axion literature. Here, we can predict the mass 
of the axion if it constitutes the main dark matter component of the 
universe. 

2. Alternatively, if Tj-eh < fa, then any strings which may have been formed 
before the inflationary epoch would be diluted by the subsequent rapid 
exponential expansion with which it is associated. In this case, we can¬ 
not, in principle, make a definite prediction about the axion mass on 
cosmological grounds. 

It is important to elaborate this ‘predictive’ distinction because, in simple 
inflationary scenarios with T^eh < fa, the estimate (|^) will still apply and a 
mass prediction appears to be meaningful. However, Oi is homogeneous over an 
entire inflationary domain, certainly exceeding the present horizon, and it is set 
to a fixed, but arbitrary, value. By (|^), this unknown value of 9, in our domain 
very precisely defines the value of rua required for the axion to be the dark 
matter—or vice versa. Either way, we have some sort of tuning re-emerging, 
which is not particularly consistent with the original motivation for the axiorj]. 
Nevertheless, in section 4 we discuss this scenario, which gives a wide range 
of possible axion masses, along with other ‘non-standard’ sceriaxios. We also 
note that topological defects may have formed during inflationtll or that they 
may have been created at the end of inflation during a phase transitior£3, 
effectively removing this new fine-tuning problem by returning us to the first 
case (i). For this reason, we describe the axion string picture as the ‘standard 
thermal scenario’; it is essentially the original axion cosmology understood 
more completely. 

If a string network does form, then its decay iuto axinns will provide the 
dominant contribution to the overall axion densityoBO, since the topological 
effects involved are much stronger than the coherent zero-momentum contri¬ 
butions of (Q). In this scenario the axion field 9 is assigned different values 
at every point in space by the topological requirement that 9 changes by 27r 
around each string. This removes the possibility of any kind of fine-tuning of 
9i, since the distribution of strings is dictated by the dynamics of the phase 
transition. However, the snec truTTi o f the radiation from these strings has been 
the subject of much debataHdQEilEj’B because it dictates the magnitude of this 

*’Linde suggests that in the infinite ‘manifold’ of a chaotic inflationary model, all values 
of 6; will be realized, soi this fine-tuning can be interpreted merely as a manifestation of 
the anthropic principle|l3; the axion-to-baryon ratio is set by 6; and humankind can only 
tolerate a narrow band of ratios! However, many find this form of ‘explanation’ decidedly 
unappealing. 
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contribution. Indeed, it appears that the over-production of axions by strings 
could almost close the allowed window of values for /a, given the lower bound 
provided by astrophysical effects. 

In recent years, we have undertaken a thorough investigation of ±he_atdiig 
radiation spectrum, using both analytic and numerical techniques 


23 


a Strir 

MI9E 


The findings were in broad agreement-with the initial work of Davis et alEH’ , 
but not with the work of Sikivie et oEjO which results in a weaker constraint. 
Taken at face value the earlier results of Davis et al. would eliminate the axion 
as a cold dark matter candidate in the standard scenario. However, although 
the underlying physics of the earlier work was correct, the model for string 
evolution on which it was based was too simplistic. Subsequent work has 
provided a more sophisticated picture of axion production by a string network 
which does leave an window open for the axion. 


1.4 Overview 

In this article, we review the cosmological constraints on the axion, in par¬ 
ticular, adding greater detail to an earlier letter on the axion string bound. 
We introduce a simple model for the evolution of a network of cosmic axion 
strings based on a ‘one-scale’ model. Using this model we calculate expres¬ 
sions for the contibution to the axion density from string loops and the long 
string network, exhibiting explicit parameter dependencies. By comparing to 
the closure density, we estimate the constraint on the Peccei-Quinn symmetry 
breaking scale /a and the axion mass Wa. In the standard scenario, we show 
that this contribution completely dominates the axion density from ‘quiescent’ 
zero-momentum production^. 

We note that there has been considerable controversy associated with the 
nature of axion string radiation and how this affects the resulting cosmological 
constraint. For this reason we provide an appendix discussing this topic in 
some detail. We present a tight mathematical argument, based on the low en¬ 
ergy Kalb-Ramond action for axion strings, which demonstrates that radiation 
goes primarily into the lowest frequency modes available. This understanding 
is shown to be in quantitative agreement with high resolution numerical simu¬ 
lations of the full U(I) field equations even at moderately high energies. This 
establishes beyond reasonable doubt the efficacy of the Kalb-Ramond action 
for describing axion string dynamics, thus providing a firm basis for the cos¬ 
mological calculations. 

Finally, we discuss ‘non-standard scenarios’ including some of the broad 
array of inflationary variants and we also consider the possibility of axion di¬ 
lution by entropy production due to the out-of-equilibrium decay of massive 
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particles. Our conclusions summarize the best current estimates of the cosmo¬ 
logical constraint on the axion. 

2 Axion string network evolution 

The evolution of axion strings is qualitatively very similar to the evolution of 
local strings due to their dynamical correspondence—as demonstrated in the 
appendix. The additional long-range field, due to the coupling to the axion, 
acts primarily to renormalize the string tension and energy density, 

M « 27r/f ln(t/(5), (2) 

where the string core width is 5 ~ f~^ and we assume the typical curvature 
radius of the strings at a time t is t. Quantitatively, on small scales i <^t, 
global strings are affected by enhanced radiation backreaction; typically in a 
cosmological context axion radiation will be three orders of magnitude stronger 
than the weak gravitational radiation produced by local strings. This difference 
will alter small-scale features such as string wiggliness and loop creation sizes, 
but not the more robust large-scale network properties. We comment further 
on the nature of string radiation in the next section. 

All the pertinent events take place in the radiation era, so for definiteness 
we consider a flat (O = 1) FRW model with 

acxt^/^ p = 32 ’ n«0.12A'r3, t« , (3) 

where p is the energy density, n is the particle number density, and Af (T) is 
the effective number of massless degrees of freedom at the temperature T. 

2.1 String network formation and the damped epoch 

The string-forming phase transition creates a tangled network permeating 
throughout the universe. The largest fraction (over 80%) of string makes 
up a random walk of long or ‘infinite’ strings, and a scale-invariant distri¬ 
bution of closed loops makes up the remainder. After formation, these strings 
experience a significant damping force due to the relatively high radiation 
background density. In the case of local strings, the dominant interaction is 
through Aharonov-Bohm type scattering. This frictional force will dominate 
the dynamics of the strings, until the Hubble damping force becomes larger. 
The temperature and time corresponding of this transition are given by 

(4) 
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where G = m~^ is Newton’s constant, fi is the string energy per unit length, 
r] is the symmetry breaking scale and t{ is the time of formation. 

There are no gauge fields present for global strings and therefore Aharonov- 
Bohm type scattering no longer pertains. In this case, one finds that Everett 
scattering is the dominant process, with a scattering cross-sectioncj 


da TT 

d9 4g[log(dg)]2 ’ 


( 5 ) 


where g ~ T is the momentum of the incident particle and S ^ f~^ is the 
width of the string. 

As for local strings, the frictional force per unit string length can be esti¬ 
mated as 

F ~ natVT^P , (6) 

where at is the transport cross-section given by 


CTt 


2g[log(g(5)]2 


T-\log{TS)]-\ 


( 7 ) 


WT ^ 1 is the thermal velocity of the particles and Ap ~ —Tv is the average 
momentum transfer per collision where T is the temperature. 

By comparison to the Hubble damping force we find that the temperature 
at which frictional damping becomes negligible is given by 


T* 

[l0g(T.//a)]^ 


Gfiiripi. 


( 8 ) 


This implies that strings oscillate relativistically and begin to radiate axions 
from the time, 




10-20 


/a \ 

IQi^GeV ) 


-4 

sec. 


(9) 


2.2 The scaling regime for string network evolution 

After the damped epoch, the strings are expected to evolve towards a scale- 
invariant regime by the formation of loops and subsequent emission of radia¬ 
tion. This scaling regime is likely to be achieved irrespective of the radiative 
mechanism, although the important parameters describing the network will be 
somewhat different. 

The overall density of strings splits neatly into two distinct parts, the long 
strings with length £> t and small closed loops with i <t, 

P = Poo + PL ■ (10) 
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The exact length scale Ic, at which the distinction between long strings and 
loops is made is unimportant since large loops are rare, but it is normally as¬ 
sumed to be comparable to the horizon. However, this distinction is important 
since loops will be radiated away quickly, whereas long strings will remain until 
domain walls form close to the QCD phase transition. 

If we ignore radiative effects, then the dynamics of the strings can be 
described by the Nambu action. The equations of motion for a string in an 
expanding background are 



--(l-X^)X, e 
a 



a 


( 11 ) 


using the temporal transverse gauge, where the string coordinates are = 
{t, X), X • X' = 0 and = X'^/(l — X^). In this gauge, the string energy is 
given by 


E = 


dae , 


( 12 ) 


where cr is a spacelike coordinate along the string. Differentiating this expres¬ 
sion with respect to time and substituting in the expression for e from (pd]), 
one can obtain 

E = -2{v^)E, p = -2{l + {v^))p, (13) 

where is the average string velocity defined by 



(14) 


with E = pV and V (x a^. In this model, {v'^) is an unknown constant to be 
determined for long strings, and {v'^) « 1/2 for string loops. Therefore, ex¬ 
cluding radiative effects and the formation of loops, the density of long strings 
and loops evolve independently according to 

Poo = -—(l + (^^))poo , PL = -—pL- (15) 

The equations (|^) do not take into account the loss of energy to loops. 
The rate of energy loss from the long string network into loops can be described 
in terms of a scale invariant production function f(£/L) where i is the loop size 
and L is the characteristic length of the long string network. This function is 
defined so that p,f{£/L)d£/L gives the energy loss into loops of size £ to £ + d£ 
per unit time per correlation volume L^. The equations for the long string 
density and the loop density are 


= --(1 + - 


—: PL = -—PL + jjf{£/L), (16) 


a 


L 



where g is a Lorentz factor and c is measure of loop production rate given by 

cp^ = ^ J\efi£/L). (17) 

In order to investigate these equations, one can substitute poo = pelt‘d 
and L = into the equation for poo- If C is constant, then we have that 

the density of long strings will scale like radiation or matter in their respective 
eras, and the evolution is self-similar that is, the large scale features of the 
string network will look the same at all times, except for a universal scaling 
proportional to the change in horizon size. Performing this substitution, one 
obtains 

^ = ( 18 ) 

where the scale factor is given by a{t) oc . This equation has a stable fixed 
point solution for C, corresponding to the scaling regime. Setting the right 
hand side of this equation to zero allows one to derive a relation between the 
loop production rate, the long string density and the average velocity of the 
long strings. 

Using a{t) oc in the radiation era, one can deduce that the loop pro¬ 
duction rate required to maintain scaling is c = — (u^)) . One can also 

solve the equation for the density of loops to give 



which can be approximated at late times by pL = p^(£t) where 

pOO 

V = / dxy/xf{x). (20) 

Jo 

Hence, the number density of loops in the interval £ to £ + d£ defined as 
li£ n{£, t) d£ = pl{£, t)d£ is given by 


n{£, t) 


V 

£bl2^Sl2 ■ 


( 21 ) 


At this stage we have to make an assumption as to the form of the loop 
production function f{x). One of the basic philosophies behind the scaling 
regime is that all the properties of the string network are constant with re¬ 
spect to the horizon. Therefore, a sensible assumption seems to be that the 
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loop production function is peaked at some constant fraction of the horizon. 
Mathematically, we treat this as a delta function 

f{x) = cS(x - , (22) 


and finally one obtains 

i. = gco^l'^CI'^ = • (23) 


These expressions do not take into account the effects of radiation on 
the strings loops which we discuss in the next section. However, this can 
be achieved simply using a linear decay of string loop energy given by ^ = 
£o — K{t — to)- Using io = at^, one can deduce that 


£ + TGgt 
l + TGg/a ■ 


(24) 


Making the substitution £ ^ /q, that is, n(£o,t)d£o ^ n{£,t)d£, the number- 
density of loops in the radiation, taking into account the decay of loop length, 
is given by 

‘ (I + 

with V redefined to be 

u = {l-\- ga^^'^C,{l — (u^)) . ( 26 ) 


All these parameters, but for the loop size_r^can be reliably estimated from 
high resolution simulations of local stringa230 (for example, ^ « 13). 


2.3 Axion mass ‘switch on’ and domain walls 

Near the QCD phase transition the axion acquires a_aiass and network evo¬ 
lution alters dramatically because domain walls formed, with each string be¬ 
coming attached to a waUPd ‘Soft’ instanton calculations give the initially 
temperature dependent masajlld 

m,(r) = (0.1±0.03)m,(^^^^j , (27) 

which achieves its maximum value for T » Aqcd at 
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This mass only becomes significant when the Compton wavelength falls inside 
the horizon, that is, m{t)t ^ 0.75 at the time 


t ~ 10.8 X 




where A is a constant of order unity which quantifies parameter uncertainties 
at the QCD phase transition |, 


A = 10 


± 0.5 


6 X 10-6eV 


0.82 


/ Aqcd 
V 200MeV 


- 0.65 


/ A/qcd N 

\ m ) 


- 0.41 


(30) 


Large field variations around the strings collapse into localized domain walls 
at t. Subsequently, these domain walls begin to dominate over the string 
dynamics when the force from the surface tension becomes comparable to the 
tension due to the typical string curvature ct ^ /r/t. 


tw ~ 1-7 X 10-®A2 


( /a m, \'Ya7qcd\ 

yW'^GeV J V6xl0-6eVy V 60 7 


0.5 


sec. 


(31) 


The demise of the hybrid string-wall network proceeds rapidljO, as demon¬ 
strated numericallyHO’Eil. The strings frequently intersect and intercommute 
with the walls, effectively ‘slicing up’ the network into small oscillating walls 
bounded by string loops. Multiple self-intersections will reduce these pieces in 
size until the strings dominate the dynamics again and decay continues through 
axion emission. 


3 The nature of axion string radiation 

Axion strings oscillate relativistically and radiate their energy primarily into 
axions; this is the preferred channel because the string coupling to the axion is 
much stronger than, say, gravitational radiation. The axion string is a global 
string with a long-range Goldstone boson field in which most of its energy 
resides, as the logarithm in (||) indicates. With /a ^ 10^° GeV this effective 
renormalization is about log(t/5) ~ 70, but the apparent ‘non-locality’ of a 
global string is more imagined than real. To illustrate this point, consider 
surrounding a straight axion string with a narrow tube of radius R ~ t/100; 
this tube would contain 95% of the string’s energy while only enclosing 0.03% of 

^Note that to calculate i, we assume an effective number of massless degrees of freedom 
JV in an epoch when its actual value is falling. However, it should be possible to use an 
averaged value throughout this epoch. 
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Figure 1: Axion radiation from an oscillating periodic string configuration in a three- 
dimensional field theory simulation. In this perpendicular cross-section, the string oscillates 
horizontally. Note the dominance of n = 2 quadrupole radiation. 


the available volume. Nevertheless, the understanding of global or axion string 
dynamics has some heuristic pitfalls and, in the past, there was considerable 
debate on the subject. In order to assure the reader of the veracity of the 
conclusions that follow on this key issue, an appendix is provided summarizing 
more recent literature on the subject. 

On cosmological scales, as we have stated already, the basic fact is that the 
axion string essentially behaves like a local cosmic string, but with stronger 
radiation damping. The radiation power arising from an oscillating string loop 
is independent of it size. This scale-invariant power is given by 

P = Ta/f = (32) 

where Fa ~ 65 has been found from numerical simulations and k essentially 
defines the radiation backreaction scale. This radiation power loss leads to the 
linear decay of the loop size with time, 

I = l\ — K{t — U). (33) 


Any loop will disappear into axions after about 10-20 oscillations. 



(a) (b) 



mode mode 


Figure 2: Power spectrum of axion radiation from an oscillating periodic string (as in fig. 2). 
Here, the initial string configuration has a sharp ‘kink’, but the initial high frequency radi¬ 
ation (a), rapidly gives way at later times to the dominant n = 2 quadrupole radiation (b). 

This spectral evolution demonstrates the effect of radiative backreaction. 

The radiation spectrum into which the loops decay is dominated by the 
lowest available frequencies, that is, the lowest harmonic wavelengths propor¬ 
tional to the loop size £. If we decompose the loop radiation spectrum into the 
power in each harmonic Pn then we expect a power law fall-off at large n, that 
is, 

P = ^P„, Pn(xn~'^ (n » 1), with q>4:/3. (34) 

n 

For typical loops, the spectral index q is expected to be greater than 4/3 
because of radiative backreaction effects. The fact that the spectrum is dom¬ 
inated by the lowest harmonics can be seen in some results from numerical 
simulations of axion strings illustrated in figs. 1-2. For this particular ax¬ 
ion string configuration the quadrupole (n = 2) is most prominent, as shown 
quantitatively in fig. 2 with a power spectrum analysis. 

The key uncertainty in cosmic axion calculations is the scale-invariant size 
£ = at at which loops are created by the string network. The parameter a is 
determined by radiative backreaction effects on the string network and, since 
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the backreaction scale is set by k in most authors anticipate a ~ k. 

However, smaller values are not necessarily excludednJ, so here we consider 
the range O.l^a /and this key parameter will appear in subsequent 
axion bounds. The above facts are all we need to know to calculate the cosmic 
axion abundance. 


4 The axion density in the standard thermal scenario 

Given the loop distribution (p5|), we can calculate the energy density of emitted 
axions. The radiation spectrum from a loop of length averaged over various 
loop configurations, is given by 

= (35) 

where the spectral function g{x) —a continuum limit of (jsj)— is normalised by 

g{x)dx = Ts., (36) 

and Ta is defined in (^. We shall assume that loops are at rest, because any 
initial velocity will be redshifted and the net error when averaged isotropically 
over all loops should be relatively small. 

The energy density of massless axions emitted at time ti in an interval dti 
with frequencies from uji to iOi + duji is 

POO 

dps^ih) = dtidivif^ / d£n{i,ti)ig{iu;). (37) 

^0 

Assuming Af constant, the spectral density can be calculated by integrat¬ 
ing over ti < t, taking into account the redshifting of both the frequency, 
Lo = a(ti)/a(t)wi, and the energy density, pa oc a“^. Neglecting the slow 
logarithmic dependence of the backreaction scale k, we have 




J a 
^ 3/2 




di 


vi 

(; + Kti)5/2 





(38) 


Under the substitution x = IjM , ^ , the range of integration is 

transformed and become^ 


^ 4/2z/ ^ ^ 


(—)) 


-3l2 


('+3 


-3/2' 


(39) 
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since the contribution from the lower limit can be shown to be zero for the 
range of oj under consideration. This implies the peak contribution to the 
density comes from those axions emitted just before wall domination. 

One can approximate the integrals of g{z) by noting that the dominant 
contribution comes in the range Att < z < 47rn*, where n* is the mode beyond 
which the radiation spectrum of loops can be truncated due to backreaction. 
Assuming 47rn* <C ujKt and a'^K, one can use the normalisation condition 
(^), to deduce that the integral (^) becomes 



4ra/f^ 


(40) 


This estimate is only formally accurate for a^K, but it should also yield useful 
information for a^K. From this expression we can obtain the spectral number 
density of axions 

_ 1 dpa (41) 

dio Lu dio 

Integrating and comparing with the entropy density of the universe, s = 
2 'k^NT^/A h, the ratio of the axion number density to the entropy at tw can 
be calculated as 



a ( K\ 3/2 

( Of \ “3/2 


6.7 X 10® 11 -h -) 

1- (l + - 

S 

V aJ 

V At/ 


- \ / -f \ 2.18 

\ / /a \ 

6xl0-6eVy \lO^^GeY) 


(42) 

(43) 


using typical parameter values Ta « 65 and v 6q!^/^( 1 + At/a)^/^. As¬ 
suming number conservation after tw and using the entropy density sq = 
2809(To/2.7K)^cm“^ and critical density pcrit = 1.88/i^ x 10“^®gcm“^ at the 
present day, one can deduce that the axion loop contribution is 


n 


a,^ 


10.7 




\ 3/2 


- 1 


h-^A 


To 

2.7K 


/a 


IQi^GeV 


1.18 


(44) 


It should be noted that the dependence of ( ^ ) on the ratio a/n comes about 
because the lifetime of a loop produced at ti is {a/K)ti. 

The contribution from long strings can also be roughly estimated assuming 
that the radiation from the long string network does not affect the scaling 
balance condition. The basis for this calculation is the radiation power per 
unit length 


dP TT^f! 
di \&'yt ’ 


(45) 
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with the long string backreaction scale given by 7 ^ ( 7 r^/ 8 )[ln(t/( 5 )]“^. As¬ 
suming the ra^ative dominance of this smallest scale "ft (as noted elsewhere 
in simulationsEj), one can calculate the spectral density of axions from long 
strings 


dpi, _ 

duj 8"fLL!t^ 


(46) 


Using similar methods to those used for loops we obtain 




1.2h-‘^A 


Tn 


2.7K 


/a 


IQi^GeV 


1.18 


(47) 


which is found to be independent of the actual backreaction scale 7 . The 
considerable uncertainty of (£7) must be emphasised given its sensitivity to 
the amplitude of small-scale structure and the overall long string radiation 
spectrum. A comparison of the two contributions (E# and (^) yields, 


n 


a,^ 




10.9 


1 + - 

K 


3/2 


- 1 


(48) 


independent of A and h. For the expected parameter range, that is 0.1 < 
aju < 1.0, the loop contribution is considerable larger. We discuss the axion 
bound quantitatively in the conclusion. 

An pxder-of-magnitude estimate of the demise of the string/domain wall 
networkEll indicates that there is an additional contribution 

(49) 


This ‘domain wall’ contribution is ultimately due to loops which are created at 
the time ^ Although the resulting loop density will be similar to (p5|), there 
is not the same accumulation from early times, so it is likely to be subdominant 
relative to (^). Both the long string and domain wall contributions will serve 
to strengthen the loop bound (^) on the axion. 


5 Alternative scenarios 

5.1 Inflation 

The discussion of the axion density presented so far presumes that a global 
string network forms after any epoch of inflation, that is, Tj-eh > /a in the 
standard thermal scenario (i) of section 1.3. In the alternative scenario when 
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Treh < fa, the axion density and any strings formed before inflation are ex¬ 
ponentially suppressed by the rapid expanslor||. In this case, the only contri¬ 
bution to the axion density comes from the initial misalignment mechanism, 
that Is, ria as given by (0). However, 9i is homogeneous on scales larger then 
the current horizon and there is no a priori reason to suppose that 9i should 
take any particular value. Essentially, with theHreedom to choose any value of 
9i, there is no constraint on /a. Some authoralj seem to favour larger values 
of 0i(^ 7’‘/2) as being more ‘natural’, since they avoid apparent fine tuning 
and anthropic arguments about our region of the universe. For example. If one 
were to chose to conservatively limit 9, to He In the range O.1^0i^7r/2, then 
the constraint on the axion becomes /a ^ lO^^-lO^^GeV , (rua ^ 0.05-50/ieV), 
which is—it has to be admitted—a rather indefinite prediction (before other 
particle physics and cosmological uncertainties are included). Further details 
about this scenario can be found in general reviewaU. 

Various attempts have been made to avoid the anthropic arguments asso¬ 
ciated with this ‘quiescent’ inflationary axion by appealing to particle physics- 
motivated models in which 9i is set by theuipnditions at the end of inflation, 
with one possibility being hybrid inflationlla. In this case, the axion held is 
coupled to the inhaton causing the formation of topological defects at the end 
of inflation, but this returns us to the standard thermal scenario described pre¬ 
viously. Such models may occur naturally In supersymmetric axion modelsu. 

5.2 Entropy production 

Entropy production by the out-of-equilibrium decay of massive particles be¬ 
tween the QCD phase transition and nucleosynthesis, can weaken all the bounds 
on /a. If the entropy is increased by some factor /3, that is s ^ /3s, then the 
axion density is decreased by a factor /3“^, that is Ha —> /3“^Ha. For example. 
It has been noted that that the dec^ of the saxino—the spin zero partner of 
the axion—can lead to a dilution byE3 

which can be up to a factor of 1000. Such entropy production appears to 
be a substantial and somewhat inelegant extension of the axion model which 

is possible to form strings during inflation^, however, the significance of this possi¬ 
bility is far from clear at the present, since quantitative predictions of the initial conditions 
are difficult to make. If just a small number of long strings survive in the initial distribu¬ 
tion, then it may be that the scaling regime can be achieved, albeit slowly. Since the most 
important axions are those emitted just before the QCD phase transition, the network will 
have plenty time to reach the appropriate scaling regime. 
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Figure 3: Parameter uncertainties in estimates of the dark matter axion mass. The key string 
parameter q/k (the loop-size/backreaction-scale ratio) is plotted from the usual assumed 
value o/k fa 1 to a possible lower limit o/k = 0.1. The bold line is the axion mass for a 
Hubble parameter h = 0.75 and the shaded region allows for errors in the range 0.5 S; h $ 0.9. 
The dotted line compounds this with particle physics uncertainties. 


combines an extra degree of fine-tuning. However, it is another mechanism by 
which to produce an axion detectable in the current search range 1-10 ^eV. 


6 Conclusions 


The cosmological axion density has been calculated in the standard thermal 
scenario by considering the dominant contribution from axion strings. Un¬ 
like the alternative scenarios described above, there is, in principle, a well- 
defined calculational method to precisely predict the mass ma of a dark mat¬ 
ter axion. For the currently favoured values of the Hubble paranaeter (Up « 
75 kms“^Mpc“^), we use string network simulation parametersE3’E®l in (44) to 
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find 


/a < 4.7 X 101° GeV > 130 ^eV, h = 0.75, (51) 

[/a < 2.3 X 10i°GeV Wa > 250 ^ieV, h = 0.5], (52) 

[/a < 6.4 X 101° GeV Wa > 78/zeV, /i = 0.9]. (53) 

The key uncertainty in this string calculation is the parameter ratio a/n (that 
is, the loop-size to backreaction scale ratio). Here, like most authors we have 
assumed a/K « 1.0, however, this remains to be firmly established quantita¬ 
tively. Gonceivably, a smaller ratio as low as a/n « 0.1 is possible and could 
weaken the bound to ma^24^eV and /a ^2.5 x lOUGeV for h = 0.75. We 
exhibit the effect of a/n on the axion mass in fig. 3. Reducing this uncertainty 
remains a key research goal, though a technically difhcult one. Also included 
in fig. 3 are the theoretical particle physics uncertainties, which are about a 
further factor of 2.5. 

Even the weakest axion string bound is considerably stronger than the 
original constraint from the ‘quiescent’ production of zero momentum axions 
for which /a ^ IQi^GeV. This bound narrows the axion window but there still 
remains a considerable range, /a ~ 10°-10^^GeV, lying above the astrophysical 
constraint. These estimates motivate the search for a slightly heavier axion 
outside the detectable range of the present axion dark matter experiment. 
Pinpointing more precisely the predicted dark matter axion mass near rua ~ 
100 /reV remains a theoretical priority. 


Acknowledgments 

EPS is particularly grateful to Professor Klapdor-Kleingrothaus and to the 
other Heidelberg workshop organisers for their invitation and hospitality. We 
would both like to thank Alex Vilenkin, Rob Caldwell, Michael Turner, Georg 
Raffelt and David Lyth for helpful discussions during this work. RAB is sup¬ 
ported by PPARC postdoctoral fellowship grant GR/K94799. This work has 
also been supported by PPARC rolling grant GR/K29272. 


7 Appendix: Axion string radiation 


7.1 The Goldstone model 


The essential features of global strings are exhibited in the simple U(l) Gold- 
stone model with action given by 


S = f d'^x 


- iA($$ - /2)2 


(54) 
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( 55 ) 


= 




where $ is a complex scalar field which has been split as $ = c/ie®’’ into a 
massive (real) Higgs component (j) and a massless (real periodic) Goldstone 
boson -d. 

For a straight global string in flat space, lying along the z-axis, the appro¬ 
priate ansatz is 

$(r, 9) = ().(r)e“^ , (56) 

where 9 is the azimuthal angle and n is the winding number, and one takes 
the usual boundary conditions, (/>(0) = 0 and (/) ^ /a as r —> oo. Despite these 
conditions, the energy per unit length ^ is logarithmically divergent. 


m(A) 



r d9 


2 'Kr dr 


Mo + 27r/f log 



(57) 


where S ^ (V^/a)”^ is the string core width and /xq /f is the core energy 
associated with the massive field </> (that is, within r^5). The length-scale A 
is the renormalisation scale provided in general by the curvature radius of the 
string or the average inter-string separation. 

It is clear from this that the string is not localised. However, in the follow¬ 
ing discussion it will be shown that the dynamics of global strings are similar 
to those of other types of strings, except for this renormalised string tension. 
Given /a 10^° — lO^^GeV, the logarithm in (|^ ) implies that there is much 
more energy in the Goldstone field than in the string core /xq- Thus the dy¬ 
namics of the string are dictated by the Goldstone self-field, a fact which has 
made the understanding of global strings intuitively hazardous. 


7.2 The Kalh-Ramond Action 

The analytic treatment of global string dynamics is hampered by the topolog¬ 
ical coupling of the string to the Goldstone boson radiation field. However, 
we can exploit the well-known duality between a massless scalar field and a 
two-index antisymmetric tensor field to replace the Goldstone boson 7} in 
( ^ ) via the relation 

= ^fae^.xpd'^B^P . (58) 

The canonical transformation generated by (|8|) requires the addition of a total 
derivative to the action. In this case the total derivative which generates the 
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transformation is given byCJ’O 


6S = 


d^xe'^‘'^Pd,B,xdpd 


d^x 


dpiB^xdpd) - Byxdpdpd 


(59) 


where is the totally antisymmetric tensor in four dimensions. Under this 
transformation the Goldstone action (^) becomes 


S = 


drx 




6(/)2 






(60) 


where the field strength tensor is = d^B'"^ + d^B^'" + d'^B^P. One 

should note that the sign of the term is the opposite of that one would 
deduce by substituting the duality relation ( p^ directly into the Goldstone 
action (|^). One can deduce the correct sign, as shown in (|^), by treating the 
antisymmetric tensor as a Lagrange multiplier. 

In the case of a spontaneously broken symmetry, the term added to the 
action (|5|) is no longer a total derivative due to the presence of a vortex at 
(/) = 0. One finds that the commutator ePP'^^dpdp'd is non zero and one can 
define 

ePP''^dpdpd = , (61) 

where is an effective current density given by 

JP’' = 1^ J ^( 4 ) _ x{a, t)) daP'^. ( 62 ) 

The area element daP'^ is given in terms of the worldsheet X {a, r) swept out 
by the zeroes of the Higgs field (</> = 0), 


daP'' = d^^daXPdbX'^dadr. 


(63) 


Rearranging the terms and evaluating the delta function in (|^) allows the 
action to be written as. 


S = 




[dp4>f 


6(/)' 







Bp.daP’'. 


(64) 


As for strings for local strings, in the neighbourhood of the string, one may 
use a local coordinate system to define. 


a;'"(cr,T,p\p^) = XP{a,T) +mPp 


a 


(65) 
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where (a, r) are the coordinates on the string worldsheet and are two 
orthonormal vectors perpendicular to the worldsheet and are related co¬ 
ordinates. If one integrates radially over the massive degrees of freedom by 
defining the bare string tension to be 


Mo = - y ^ ^ 

one can deduce the Kalb-Ramond action@ 




S = -fio J y/^dadr + ^ J H'^ - 27:fa J B^yda^’'. 


( 66 ) 


(67) 


where ^ab = gfiudaX^dbX'^ is the induced metric on the worldsheet and 7 = 
det( 7 ab). The first term is the well known Nambu action of local strings, the 
second is the antisymmetric tensor field strength for both external fields and 
the self field of the string and the last is a contact interaction between the 
field and the string worldsheet. This coupling between the string and 
Bfj^i, is analogous to the electromagnetic coupling of a charged particle, and 
is amenable to the same calculational techniques. This is the basis for the 
subsequent analytic work reviewed here. 


7.3 Equations of motion 

Varying the Kalb-Ramond action (^) with respect to the worldsheet coor¬ 
dinates and the antisymmetric tensor, gives the equations of motion for the 
string and the antisymmetric tensor field equation, 

goda iV^r’^dbXf^) =E>^ = , (68) 

= _47rJ“/3 = -27: fa J dadTS^{x - V(ct,t)) e’^^daX^dbX^ . (69) 

One can use the conformal string gauge -|- X'^ = 0 and X ■ X' = 0, and 
also the Lorentz gauge of the antisymmetric tensor field, that is, dfj,B^‘^ = 0, 
to deduce that 

- X"^) =E>^ = 27:faH^^‘^d(^Xa,x'f^ - X'a^Xp) , (70) 

UBa ,0 = 27: fa J dadriXo^X'p - V^,X^)5(4) _ x{a, r)) . (71) 

As for the point electron one can split up the force into two parts; one due 
to the self-field and the finite radiation backreaction force. The equations of 
motion are therefore 

Aio (V'" - X"'^) = , (72) 
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Figure 4: Loop and periodic long string trajectories which have been the subject of extensive 

analytic and numerical studies. 


where the self-force is 


^r = -27r/flog(A/<5) 


X>^ - 


(73) 


and a fipt-nrHpr approximation to the hnite radiation backreaction force is 
given bycilld 


JPf" = 


4^/fA 



X>^ + 


(x'-x 

X2 


x^' 


(74) 


where A is the renormalisation scale introduced in (|^) . The self field is loga¬ 
rithmically divergent, but it can be absorbed such that the equations of motion 
in the conformal gauge becomeQ^ 

^(A) {X^ - A"^) = . (75) 

The renormalised equations of motion (fz^ ) for the string can be approximated 
by the Nambu equations of motion, assuming the effects of radiation backre¬ 
action to be small, that is « 0. 

These equations have closed loop and periodic long (or infinite) string 
solutions. The loop solutions are parametrized by their invariant length of 
the loop L, which is closely related to the characteristic frequency 17 = 27r/L, 
whereas the long periodic solutions are parametrized by their wavelength L and 
the ratio of amplitude to wavelength or the relative amplitude £ = 2 ttA/L^ 
where A is their amplitude. Fig. ^ shows a schematic of the two types of 
solution we have considered in detail. 
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7.^ Appendix conclusion: The nature of string radiation 

The nature of radi^aaTLom oscillating axion (and local) strings has been ex¬ 
tensively studiedOEJ’OLd. The methods developed allow one, in principle, to 
calculate the radiation power from any arbitrary periodic loop or long string 
configurations, yielding the magnitude of the total power as well as the spec¬ 
trum of the radiation. For arbitrary long string configurations we have also 
derived a linearized expression whictL^ows one to calculate the leading order 
contribution to the radiation poweiollj. By thorough examination of the ex¬ 
pressions for the radiation power and applying them to various loop and long 
string configurations, we have been able to make the following deductions: 

(i) The radiation power for loops is independent of the loop sizeSEl with the 
power given by P = Ba/a = K/r which leads to the linear decay in the loop size 
i = Iq — n{t — to) as in (^. 

(ii) Loops radiate principally in their fundamental mode with a radiation spec¬ 
trum such that oc n~‘> for large n where q > 4/3. 

(iii) The radiation from long strings is generically in the fundamental mode. 
However, for exactly periodic solutions the fundamental mode is suppressed 
and the radiation is quadrupole (as in figs. 1-2). 

(iv) The radiation power from long string configurations is proportional to , 
which leads to the exponential decay of their amplitude £. Again, for exactly 
periodic configurations cancellation gives radiation proportional to E'^, leading 
to a power-law type decay. 

(v) Some string configurations, like kink solutions, which theoretically have a 
pathological divergence in their radiation power (P„ oc n~^), rapidly revert to 
a smooth spectrum dominated by the lowest harmonics because of radiation 
backreaction effects (see fig. 2). 

These predictions broadj^mphold the assumptions behind the original work on 
global string radiationEElBB but Idiai axe completely contrary to the predictions 
of a ‘harder’ Pn oc n~^ spectrumEdO’Ed. We conclude that axion radiation from 
axion strings, when extrapolated to cosmological scales, will be very similar to 
gravitational radiation from Nambu strings. 
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